Inequivalence of ensembles in a system with long range interactions 
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We study the global phase diagram of the infinite range Blume- Emery-Griffiths model both in the 
canonical and in the microcanonical ensembles. The canonical phase diagram is known to exhibit 
first order and continuous transition lines separated by a tricritical point. We find that below the 
tricritical point, when the canonical transition is first order, the phase diagrams of the two ensembles 
disagree. In this region the microcanonical ensemble exhibits energy ranges with negative specific 
heat and temperature jumps at transition energies. These results can be extended to weakly decay- 
ing nonintegrable interactions. 
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Systems in d dimensions with a pairwise interaction po- 
tential which decays at large distances as V{r) ^ i/j-d+a 
with —d < cr < 0, are referred to as nonintegrable, or 
systems with long range interactions. Such systems have 
an ill defined thermodynamic limit [Q. This may be 
correctly restored by applying the Kac prescription 
within which the potential is rescaled by an appropriate, 
volume dependent, factor which vanishes in the thermo- 
dynamic limit. However, even within this scheme, the 
energy remains nan additive, i.e. the system cannot be 
divided into independent macroscopic parts, as is usually 
the case for short range interactions. This fact has no 
dramatic consequences if one is restricted to the canoni- 
cal ensemble, but it produces striking phenomena in the 
microcanonical ensemble. For example it may result in 
a negative specific heat, as was first clearly discussed by 
Lynden-Bell [| and Thir ring MM- Indeed, it has been 
originally observed by Antonov ^ that classical gravi- 
tational systems (a = —2, d = 3) show features of such 
kind. However, here the physical situation is made more 
complex by the presence of a singularity of the interac- 
tion potential at short distances. For a careful discussion 
of the statistical mechanics of these systems see Ref . [D . 

In the present Letter we consider a simple model for 
which the main features of the phase diagram can be 
derived analytically both within the canonical and the 
microcanonical ensembles. We demonstrate that in the 
region where the phase transition in the canonical ensem- 
ble is first order, the two ensembles are not equivalent, 
yielding two distinct phase diagrams. The model we con- 
sider is the Blume-Emery-Griffiths (BEG) model with in- 
finite range interactions (cr = --d). This is the simplest 
model known to exhibit both continuous and first order 
transition lines. It is defined on a lattice (hence, diver- 
gences at short range are removed), where each lattice 
point i is occupied by a spin-1 variable Si = 1,-1,0. 
The Hamiltonian is given by 



N T ^ 

i=l 1=1 

where J > is a ferromagnetic coupling constant and 
A > controls the energy difference between the mag- 
netic {S = ±1) and the non-magnetic {S = 0) states. 
Each spin interacts with every other spin and the cou- 
pling constant J is scaled by to ensure extensivity 
of the energy. This is just the Kac prescription applied 
to our model. However, this does not entail additivity, 
in the sense that for a system made of two parts, X and 
Y, such that Hx+y = Hx + Hy + Hxy: the Hxv inter- 
action term never becomes negligible in the thermody- 
namic limit. This property applies to all thermodynamic 
potentials. 

The canonical phase diagram of this model has been 
studied in the past [^. At T = the model exhibits a 
ferromagnetic phase for 2 A/ J < 1 and a non-magnetic 
phase otherwise. The {T, A) phase diagram displays 
a transition line separating the low temperature ferro- 
magnetic phase from the high temperature paramagnetic 
phase (see Fig.0). 

The transition is first order at high A values and be- 
comes continuous at low A. The critical (second order) 
line is given by 

pj ^ ^e^^ + 1 , (2) 

where f3 = l/ksT. The two segments of the transition 
line (high and low A) are separated by a tricritical point 
located at A/J = ln(4)/3 ~ 0.4621, f3J = 3. The first 
order segment of the transition line is obtained numeri- 
cally by equating the free energies of the ferromagnetic 
and the paramagnetic states. 

We now consider the phase diagram of the BEG model 
(|l|) within the microcanonical ensemble. Let N^, iV_, A^o 
be the number of up, down and zero spin, respectively, in 
a given microscopic configuration. Clearly, A^+ -I- N- + 



1 



A^o = N. The energy E oi a configuration is obviously a 
function only of N+, N- and Nq. It does not depend on 
the specific spatial distribution of the spin variables. It 
is given by 



E = AQ- 



_J_ 

2N' 



(3) 



where Q = J2iLi ~ ^+ + quadrupole mo- 

ment and M — — ~ magnetization 

of the configuration. In order to calculate the entropy of 
a state with energy E we note that the number of mi- 
croscopic configurations compatible with macroscopic 
occupation numbers , 7V_ and A''o is 



m 



N+\N-\No\ 



(4) 



Thus, in the large N limit, the entropy 5 = /c^ Inil cor- 
responding to these occupation numbers is given by 



S = -kBN[{l - q) ln(l - q) + ^{q + m) \n{q + to) 



-((7 - m)ln(g - m) - 9ln2] , 



(5) 



where q — Q/N and to = M/N are the quadrupole mo- 
ment and magnetization per site, respectively. 

Let e — E/AN be the energy per site, normalized by 
A. Equation (|) may be written as 



q = e + Km , 



(6) 



where K = J/2 A. Using this relation, one expresses the 
entropy per site s = S/{kBN) as a function of to and e. 
Maximizing s{e,m) with respect to m one obtains both 
the spontaneous magnetization ms(e) and the entropy 
s(e) of the system for given energy. In order to locate the 
continuous transition line between the paramagnetic and 
the ferromagnetic phases we expand s(e, m) in powers of 
TO. This expansion takes the form 



.s = So + Arn'^ + Brn^ + 0{rri^) , 



(7) 



where sq = ^(e, to = 0) is the entropy at zero magnetiza- 
tion 

so = -(l-e)ln(l-e)-elne + eln2 , (8) 
and A and B are the expansion coefficients 
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In the paramagnetic phase both A and B are negative, 
and the entropy is maximized by m = 0. The continuous 
transition to the ferromagnetic phase takes place sX A = 



for _B < 0. In order to obtain the critical line in the 
(T, A) plane we note that the energy e is related to the 
temperature by the usual thermodynamic relation 



A 



ds 



knT de 
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Making use of the fact that the magnetization to vanishes 
on the critical line one obtains 



A 



In 



2(1 -e) 



(11) 



This relation, together with the equation A — 0, yields 
the following expression for the critical line 



2PK =^e^ + 1 



(12) 



where f3 = /3A. Equivalently, this expression may be 
written as PK = l/2e. The microcanonical critical line 
thus coincides with the critical line (||) obtained for the 
canonical ensemble. The tricritical point of the micro- 
canonical ensemble is obtained at A = B = 0. Com- 
bining these equations with ( pi] ) one finds that at the 
tricritical point P satisfies 
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Equations (|1^ Jl^) yield a tricritical point at X ~ 1.0813, 
(3 ~ 1.3998. This has to be compared with the canon- 
ical tricritical point located a.t K = 3/ln(16) ~ 1.0820, 
(3 = ln(4) ~ 1.3995. It is evident that the two points, al- 
though very close to each other, do not coincide and the 
microcanonical critical line extends beyond the canoni- 
cal one. In the region between the two tricritical points, 
the canonical ensemble yields a first order transition at a 
higher temperature, while in the microcanonical ensem- 
ble the transition is continuous. 

To study the microcanonical phase diagram we con- 
sider the temperature-energy relation Tie). This curve 
has two branches: a high energy branch (0) correspond- 
ing to TO — 0, and a low energy branch obtained from 
( pO| ) using the spontaneous magnetization TOs(e). At the 
intersection point of the two branches the two entropies 
become equal. In Fig. |^ we display the T(e) curve for 
increasing values of A. For A/J = ln(4)/3), correspond- 
ing to the canonical tricritical point, the lower branch 
of the curve has a zero slope at the intersection point 
(Fig. ^). Thus, the specific heat of the ordered phase 
diverges at this point. Increasing A to the region be- 
tween the two tricritical points a negative specific heat 
in the microcanonical ensemble first arises {dT/de < 0), 
see Fig. |^. At the microcanonical tricritical point A the 
derivative dT/de of the lower branch diverges at the tran- 
sition point, yielding a vanishing specific heat (Fig. ^). 
For larger values of A a jump in the temperature appears 
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at the transition energy (Fig. ||d). The lower tempera- 
ture corresponds to the m — sohition ( |ll|) and the up- 
per one is given by exp(/3) — 2(1 — q*)/\/ (q*)^ — {m*Y, 
where to* , q* are the values of the order parameters of 
the ferromagnetic state at the transition energy. The 
negative specific heat branch disappears at larger values 
of A, leaving just a temperature jump (see Fig. ^). In 
the A/J — !■ 1/2 limit the low temperature branch, corre- 
sponding to g = m = 1 in the limit, shrinks to zero and 
the TO = branch (|l^) describes the full energy range 
(Fig. ^). In the inset of Fig. |l| we report the transi- 
tion temperatures in the microcanical ensemble against 
A/J for both the m — Q (lower dot-dashed line) and the 
TO / solutions (upper dot-dashed fine). The fines are 
drawn starting at the canonical tricritical point. The re- 
gion between the two tricritical points is too small to be 
appreciated in the figure. A schematic phase diagram in 
the first order region is given in Fig. ^, where we ficti- 
ciously expand the region of the tricritical points. Note 
that the canonical first order line necessarily crosses the 
upper microcanonical transition line at some point. 

That such unusual effects in the microcanonical en- 
semble are associated with a first order canonical phase 
transition was also suggested in Ref . ||l^ . These authors 
discuss also short range interactions, for which such fea- 
tures are produced by finite size effects. 

As usual for mean- field models, one can express the 
free energy /(T, to) in the canonical ensemble as a func- 
tion of T and m. The spontaneous magnetization ms(T), 
the temperature-energy relation T{e) and the free en- 
ergy /(T) may be obtained by minimizing f{T,m) with 
respect to to and using well known thermodynamic rela- 
tions. We now note that the negative specific heat branch 
of the microcanonical ensemble corresponds to a local 
maximum of the free energy f(T,m) with respect to to. 

This result can indeed be derived on quite a general 
ground. It is easy to show that an extremum of f(T, to) 
corresponds to an extremum of s{e,m) with respect to 
TO. Indeed, the free energy f(T,m) may be obtained 
by minimizing fiT^e^m) = Ae — s/[3 with respect to e, 
keeping T and to fixed. This minimization yields the 
usual temperature-entropy relation ([l0|). Further mini- 
mizing /(T, e, to) with respect to to yields the result that 
9/(T, to)/9to and ds{t,m) / dm are proportional to each 
other and thus vanish together. It can be shown |l|] by 
studying the second derivatives that when the stationary 
point of /(T, e, to) with respect to e and to is a saddle 
point, the resulting entropy exhibits a negative specific 
heat. As a consequence, we can recover the full micro- 
canonical solution by studying the stationary points of 
the function /(T, e, to). However, this function is not 
typically available for non mean-field models. 

The relevant features of the BEG model with infinite 
range couplings persist also for nonintegrahle interac- 
tions. In order to investigate this point, we introduce 
a generalization of the BEG model given by the Hamil- 



tonian 

H = ^ts'. - jT.'-§ ^ (14) 

i—l i>j ^3 

where is the distance on a ID lattice between spins 
at sites i and j. The interactions are non integrable for 
a < 1. The normalization TV = 2"(7Vi-" - 1)/(1 - a) 
ensures that the energy is extensive. Models of this kind 
have been previously introduced by other authors, and 
studied within the canonical ensemble . We apply pe- 
riodic boundary conditions (p.b.c), for which the model 
is more easily tractable, and we then take r^j to be the 
smallest of the two distances compatible with p.b.c. The 
interaction matrix {rij)^°'/N can be exactly diagonal- 
ized, which allows to solve model (|lj) in the canonical 
ensemble. When appropriately rescalcd thermodynamic 
quantities are chosen, the solution is the same as for 
the a = case (as it happens for the models studied 
in Ref. fl^). Moreover, using a Fourier representation 
of Hamiltonian ( |l^ ) and considering only the long wave- 
length components, it is possible to obtain an approx- 
imate expression for the entropy in the microcanonical 
ensemble ||ll|]. Maximizing this expression at fixed en- 
ergy, we find that the a = microcanonical solution is 
also left unchanged. Since both the canonical and the mi- 
crocanonical solutions are not modified, we conclude that 
ensemble inequivalence persists for the slowly decreasing 
case a < 1. Details of this analysis will be reported else- 
where [ pl] |. 

In summary, we have compared the canonical with 
the microcanonical solutions of the infinite range Blume- 
Emery-Griffiths model. We find that the global phase 
diagrams are different in the two ensembles. Although 
they are found to be the same in the domain where the 
canonical transition is continuous, they differ from each 
other when the canonical transition is first order. Nega- 
tive specific heat and temperature jumps at the transition 
energy are found in the microcanonical ensemble. These 
results generalize those of Ref. Q in the context of a 
simple model, where by varying a single parameter one 
can observe a variety of possible features of the phase 
diagram. Moreover, we are able to understand the role 
played by the constraint of fixing the energy in the mi- 
crocanonical ensemble, which produces a stabilization of 
canonically unstable solutions. In the phase coexistence 
region, the unusual microcanonical thermodynamic prop- 
erties should result in some peculiar dynamical behavior, 
as has been observed in studies of a different mean-field 
model with continuous variables jl^ . Our results for the 
BEG model are not limited to the infinite range case, 
but can be extended to weakly decaying nonintegrahle 
interactions. 
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FIG. 1. Transition lines in the canonical ensemble. The 
critical line (dotted) ends at the tricritical point (•), where 
the transition becomes first order (full) . The first order region 
is zoomed in the inset, where we show again the canonical 
first order line (full) and the microcanonical transition lines 
(dot-dashed) 
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FIG. 2. Temperature versus energy relation in the mi- 
crocanonical ensemble for different values of A. The dotted 
horizontal line in some of the plots is the Maxwell construc- 
tion in the canonical ensemble and identifies the canonical 
first order transition temperature. 




FIG. 3. A schematic representation of the phase diagram, 
where we expand the region around the canonical (CTP) and 
the microcanonical (MTP) tricritical points. The second or- 
der line, common to both ensembles, is dotted, the first order 
canonical transition line is full and the microcanonical transi- 
tion lines are dashed (with the bold dashed line representing 
a continuous transition). 



4 



